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ABSTRACT 

A new predictive control approach is employed to control a non-linear 4 
degrees of freedom vehicle suspension model including pitch and bounce 
motions. The front and rear suspension actuator forces are considered as 
control inputs. These are determined by minimizing a performance index 
defined as a weighted combination of conflicting objectives, namely ride 
comfort and road holding. In this way, the suspension deflection and the 
control forces should be restricted to admissible practical ranges. In this paper, 
an optimization process is used to develop a multi-input multi-output (MIMO) 
non-linear control law in a closed-form which is mathematically tractable and 
suitable for implementation. To investigate special features of the derived 
control law, the closed-loop system dynamics are first evaluated analytically 
and then the corresponding simulation results are presented. The results show 
that a compromise between mentioned objectives and control energy, even in 
the presence of uncertainties, can be easily made by regulating the control 
weighting factors. 

Keywords: Non-linear suspension system; predictive control; pitch and bounce 
motions; optimization; MIMO system 


I. INTRODUCTION 

Vehicle suspension systems play an important role in the performance of modern 
road vehicles. They are mainly designed to isolate the vehicle and its occupants 
from the road roughness. This objective is to comply with the ride comfort 
criterion. On the other hand, a suspension system must ensure good contact 
between the tires and the road to meet the road holding criterion [1,2], The road 
holding is important because it ensures vehicle safety on rough roads during 
turning, braking and accelerating. The two mentioned objectives are in conflict 
with each other and a tradeoff between them is necessary in the design of a 
suspension system. Another design factor is the restriction of suspension 
deflection. Large suspension deflection causes excessive pitch and roll motions in 
abrupt manoeuvers. 

A passive suspension consists of constant elasto-damping elements and only 
provides uniform behaviour. Therefore, it cannot efficiently compromise between 
the ride comfort and handling in different situations [3], Designing an active 
suspension system despite its cost and complexity can highly compensate for the 
drawbacks of passive systems. One of the main constraints in the design of an active 
system is power limitation of its actuator, which is required to apply large forces in 
short durations [4, 5], This constraint should be considered as an important design 
factor in the active suspension system [6]. 


* Corresponding author e-mail: mirzaei@sut.ac.ir (M Mirzaei) 


Vol.34 No. I 2015 



Non-Linear Predictive Control of Multi-Input Multi-Output Vehicle Suspension System 


The first step in the design process of an efficient active suspension system is to 
select a suitable model referred as design model. The design model should include 
the main features of a real car suspension system. Various suspension system 
models in terms of number of degrees of freedom (DOF), complexity and design 
objective have been used by different researchers. Among them, the quarter-car 
model is the most common-used design model due to its simplicity. However, it only 
captures the vertical motion of suspension and does not include the pitch motion 
which becomes very important in vehicle manoeuvers. Therefore, using a half-car 
suspension model considering the pitch motion is preferred as an appropriate model 
in studying the behaviour of road vehicles. It should be noted that, designing the 
controller based on a comprehensive vehicle model with large degrees of freedom 
is often mathematically complicated. 

Another important factor to achieve a more accurate model is choosing realistic 
elasto-damping models for suspension elements. Kim and Ro [7] compared the 
responses of linear and non-linear suspension models for different road inputs and 
concluded that the non-linear characteristics of suspension elements should be 
considered in developing a more accurate system model. There are many studies 
that have used a half-car model but with linear elements [8-11]. Generally, a non¬ 
linear half-car model considering pitch and bounce motions needs to be developed 
for the design of active suspension systems. 

Beside the dynamic models highlighted above, different control strategies have 
been reported for active suspension design. Commonly used control approaches of 
suspension systems include sliding mode control [7, 12-14], non-linear adaptive 
control [9, 15], fuzzy logic control [8, 15, 16, 17], back-stepping control [18, 19], 
multi-objective robust control [20, 21], H v control [4, 22, 23], LQ-based control 
[24], LQG control approach [25, 26] and feedback linearization control [27], Each 
of these control methods considers some aspect of active suspension problem. 
However, as mentioned before, the main control problem for the active suspension 
involves finding a suitable external force that trades off between conflicting 
objectives. Additionally, the control law must minimize the required injected forces. 
Thus, an optimal control problem should be considered to find the control law by 
minimizing the weighted combination of conflicting objectives as performance 
index. On the other hand, the use of a non-linear model for suspension elements can 
broaden the valid range of operation as well as improve the accuracy of the model. 

Classic optimal control strategies based on linear models have been widely used 
to design active suspension controllers [6, 10, 28-30], Although these methods are 
systematic and can provide a good compromise between conflicting objectives and 
control force, they are exclusively applicable to linear models, whereas a real 
vehicle suspension has intrinsic non-linear behaviour. There are some non-linear 
control approaches which have used the optimization in controller design [31-33], 
However, they have used numerical or search methods for calculating the control 
input at each instant of control process. In general, the application of classic optimal 
control theories for non-linear systems requires solution of the derived non-linear 
differential equations. There are no analytical solutions for all these equations. In 
addition, numerical computation methods require online dynamic optimization 
which is difficult to implement [34, 35]. 

In this study, a predictive control approach is employed to design a non-linear 
MIMO controller for active suspension system using an optimization process. The 
method directly employs non-linear elasto-damping elements and analytically leads 
to a closed form control law which is mathematically tractable and suitable for 
implementation. In addition, online numerical computations are not required at each 
sample time. Following the success of this method fora quarter-car suspension 
model which is related to a single input system [2], the authors decided to extend 
this approach for a realistic model considering pitch and bounce motions. In this 
case, a MIMO control problem with two inputs is solved analytically. To investigate 
the special features of the derived control laws, the closed-loop system dynamics are 
evaluated and simulation results are presented. 
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Figure 1.4D0F half-car model and its free body diagram. 


2. SUSPENSION MODEL 

In this section, a 4DOF half-car model containing essential features of the vehicle 
suspension system is proposed for designing the controller. The bounce and pitch 
motions of the sprung mass together with the front and rear wheels hop constitute 
the DOF for this model (Figure 1). The equations of motion for the model are simply 
derived using Newton’s second law as follows: 

For the sprung mass: 


m s x s =-ff-f r+ u f+ u r (1) 

Id = -af f + bf r + au f -bu r (2) 

and for the un-sprung mass 

m f x f = ff~ fjt~ u f (3) 


m r x r =f r -fn-u r (4) 

where x s is the sprung mass displacement, x f and x r are the front and rear tire 
displacements, respectively, and f r represent the front and rear elasto-damping 
forces, respectively. Also, m s , ra^and m r denote the sprung mass, the front un-sprung 
mass and the rear un-sprung mass, respectively.^ and/ n are the visco-elastic forces 
of front and rear tires, respectively. u f and u r represent the front and rear actuator 
forces, respectively, a and b are respectively the front and rear suspension distances 
from vehicle center of mass, I is the sprung mass moment of inertia about the body 
centroidal axis and 6 is the pitch rotation angle of the sprung mass. 

In linear suspension models, the spring and damper forces are often described as 
linear functions of strut deflection and relative velocity. However, the linear 
behaviour of suspension spring and damper is valid only for small displacements 
from their equilibrium position. For example, the spring has a hardening behaviour 
for deflections beyond the linear region. Kim and Ro [7] proposed polynomial 
functions for elasto-damping forces. They obtained these functions by fitting 
polynomial curves as follows to the measured data: 

/„ = k x Ax+k 2 Ax 2 + k 3 Ax 3 (5) 

f d = Cl Ax + c 2 Ax 2 (6) 

where and f d are spring and damper forces respectively, and Ax and Ax represent 
suspension deflection and velocity, respectively. k ] , k 2 and k 3 are stiffness 
coefficients and c x and c 2 are damping coefficients. 
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The state variables are defined 


x l = x s +ad - x f 

(7) 

x 2 =x s +ad 

(8) 

x 3 =x s -bd-x r 

(9) 

x 4 =x s -bO 

(10) 

x 5 =x f - Xrf 

(ID 

X 6 = X r~ X ,rr 

(12) 

X 1 = x f 

(13) 


x s =x r (14) 

where x x and x 2 are the suspension deflection and absolute velocity of sprung mass 
in front station, x 3 and x 4 are the suspension deflection and absolute velocity of 
sprung mass in rear station, x 5 and x 6 are the front and rear tire deflections, and x 7 
and x 8 are the front and rear tire velocities, respectively, and x rr are the road 
elevation in front and rear stations, respectively. 

Now, by defining the above state variables, the front and rear suspension and tire 
forces are obtained as: 


f f =k f3 xl+k f: x 2 l +k fl x l 

+ c f2 (x 2 -x 7 ) + c fl (x 2 -x 7 ) 

(15) 

f r =k r3 x l+k r2 xl+k rl x 3 

+ c r2 (x 4 -x 8 ) 2 + c rl (x 4 -x s ) 

(16) 

f ft = k “f X 5 

+ C “f( X 7 ~ X <f] 

(17) 

fn = KrX ( , 


(18) 


where , k^ 2 , k f3 and , cy 2 are stiffness and damping coefficients for front elasto- 
damping elements, k rl , k r2 , k r3 and c rl ,c r2 are stiffness and damping coefficients for 
rear elasto-damping elements. The coefficients k^ and c uf denote the stiffness 
and damping of the front tire and the coefficients k ur and c ur denote the stiffness and 
damping of the rear tire, respectively. 

The road input for the rear wheel is a delayed form of the road input for the front 
wheel: 

Xrr-X^t + X) (19) 

Here, r is the time delay between the front and rear wheels in crossing over the 
road and depends on the wheel base and speed of the vehicle: 
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T = — (20) 

where v is the vehicle speed. Now, the state space equations are obtained by 
substituting the state variables from eqns. (7) - (14) into eqns. (1) - (4): 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

where the terms g- (i = 2,4,7,8) and d- (i = 1, ..., 4) have been defined as: 


(-f f -f r ) + j(-f f a + f r b) 

(29) 

(-f f -f t )-j(~f f a + f r b) 

(30) 

m f 

(31) 


(32) 

, 1 a 2 

d.= — + — 
m s I 

(33) 

d 2=d3 = ±_± 

2 3 m s I 

(34) 

^ 4 =- + - 
m s I 

(35) 
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3. CONTROL SYSTEM DESIGN 

In active suspension systems, the controller is designed to make a suitable trade-off 
between the conflicting objectives with the restricted control forces. In this study, the 
conflicting objectives are the front and rear sprung mass accelerations, x 2 and x 4 , 
the front and rear suspension deflections, x 1 and x 3 , and the front and rear tire 
deflections x 5 and x (> . The mentioned trade-off between these objectives could be 
achieved by employing an adjustable control law which is presented and discussed in 
the next sections. 

3.1. Development of the control law 

In this section, an optimal control law is derived based on predicting the suspension 
system responses from the non-linear state equations of the half-car model. 
According to the suspension system requirements described previously, the system 
outputs are assumed as: 



(36) 


It should be noted that the sprang mass absolute velocities, x 2 and x 4 are 
considered as controlled variables and a measure of ride comfort instead of the 
sprang mass acceleration [2, 6, 7, 26, 29, 30], This is because the system degree 
relative to the sprang mass acceleration is zero, whereas, the relative degree of each 
predicted output has to be at least one. It will be shown in simulation studies that the 
sprang mass acceleration is perfectly controlled by controlling the absolute velocity. 

Now, according to the proposed control method, the non-linear response of each 
output, x. (i=l,...,6), for the next time interval, x x (t + h), where h denotes the 
prediction time, is first predicted by a Taylor series expansion and then the current 
control forces, u(t) = [uy w r ] r , are found based on continuous minimization of the 
predicted responses. 

To achieve the compromise between conflicting objectives, a performance index 
which penalizes the regulation error of output responses in the next instance and the 
current control expenditure is defined as: 



(37) 


where p ; > 0 and p 7 and p g a 0 are weighting factors indicating the relative 
importance of the corresponding terms. The predicted response for each output in 
the next interval, xft + h) (i=l, ... ,6), is approximated by a k ; th-order Taylor series 
expansion at t as follows: 




(38) 


The expansion order, k ( , for each output is selected to be suitable to design the 
controller based on prediction. To achieve a small control effort and to prevent 
complexity in derivation and implementation of the control law, the expansion order 
is limited to the relative degree of the non-linear system [34, 35], This selection 
which is related to zero control order, eliminates the derivatives of the control input 
in the prediction of each output and results in adequate performance for the non¬ 
linear system with a lower relative degree [2, 13]. Thus, the control forces will be 
constant in the prediction interval: 


— u f (t + r) = 0 /or r g[0,/i] 


( 39 ) 
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~^u r (t +t) = 0 for r g[0,/z] (40) 

The relative degree is determined as the lowest order of the derivative of each 
output, xft) (i = 1, in which the input u^(t) or lift) appears explicitly. 

According to the state space eqns. (21)-(28), the relative degree of freedom with 
respect to x 1 , x 3 , x 5 and x 6 is r = 2 and with respect to x 2 and x 4 is r = 1 . Thus, a first- 
order Taylor series expansions are sufficient for x 2 and x 4 but a second-order Taylor 
series expansions are required for x 1 , x 3 , x 5 and x (y So, the non-linear response of 
each output for the next time interval is predicted as follows: 

x i (t + h)=x i (t) + hx i (0 + ^(0 (z = 1,3,5, 6) (41) 

x i {t + h)=x i (t) + hx i {t) [i = 2,4) (42) 


Substituting the state space equations and state variables derivatives into the 
corresponding xft + h ) yields: 

{t + h) = G ; + l t u f + m,.«, (i = 1,.. .,6) (43) 

where Gs, s and ms (i = 1, ...,6) are as follows: 


+ h[x 2 -x 7 )+ t ^(g 2 -g 7 ) 

(44) 

G 2 =x 2 + hg 2 

(45) 

f+ h(x 4 -x s ) + *t(g 4 -g s ) 

(46) 

G 4 = x A + hg A 

(47) 


(48) 

+ h(x i -x rr ) + ^(g s -x rr ) 

(49) 


(50) 

h 2 d 2 

m. = -- 

1 2 

(51) 

/ 2 = hd x 

(52) 

m 2 = hd 2 

(53) 
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l - k % 

(54) 

3 2 " 


h 2 ( , ^ n 

(55) 

/ 4 = m, 

(56) 

m 4 = M 4 

(57) 

h = ~ — 

(58) 

2 m f 


m 5 = 0 

(59) 

4=0 

(60) 

h 2 

m, = - 

(61) 

2m r 


After substituting eqn. (43) into the performance index (37), the following 
necessary conditions for optimality are applied: 

£—0 (62) 

du f 

"-0 (63) 

du r 


2 A ( G * + l i u f + m Mr)h + Pi u f = 0 

2 A ( G , + h u f + + p s u r = 0 

Solving the above equations for t^ and u r yields: 

u e^PiG^-e^PiG-X 


|PjGj rn i 


(64) 

(65) 


( 66 ) 


( 67 ) 
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where e v e 2 and e 3 are defined as follows: 

=%Pi l i + Pi 


( 68 ) 


(69) 




3.2. Evaluation of the control laws 

The main contribution of the present study is derivation of the analytical optimal 
solution for the MIMO control problem of non-linear active suspension system. In 
this section, the main features of the derived closed-form control laws (66) and 
(67) which are mathematically tractable and suitable for implementation are 
described further. 

As mentioned before, the control purpose is to compromise between some 
conflicting objectives. The importance of each objective is emphasized by assigning 
a higher amount of its weighting factor relative to the others. Special cases of the 
control laws (66) and (67) can be obtained when only a specific objective is to be 
taken into consideration. Therefore, by setting other weighting factors to zero except 
p i and then by substituting eqns. (68)- (70) into the general case of the control laws 
(66) and (67), the following especial versions of the control laws are derived. 


PsPih 

U, - / \ G, 

Pi{p,m, +p s )+p s p i l i 

{i = h. 

..,6) 

(71) 

u r =- 7 2 G ; 

P,[pJi + Pt ) + Pi P,m t 

(i- 1 

,..,6) 

(72) 

To investigate more specific forms of eqns. (71)-(72), the ride comfort criterion 
i.e. minimizing the responses xfi = 2,4), as a case study are considered in details. 
However, other responses could be analysed in a similar way. Substituting 
corresponding G.s from eqns. (44)-(49) into eqns. (71) and (72) yields: 

u f=-ji x i +h Si) 

O'= 2,4) 


(73) 

U r=-j{ X ,+ h 8.) 

(i- 2,4) 


(74) 


where /c and \ are definedby substituting corresponding l t s and m i s from eqns. (50)- 
(61) into eqs (71) and (72): 


Pi{Pih 2 df +p s ) + p g p i h 2 df_ l 

_ _ PiP,h 2 f, _ 

Ps[p i h 2 dl l + p,) + p.p^df 


- (i = 2,4) 

(75) 

a-2,4) 

(76) 
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The special version of the control laws (73) and (74) are based on only the ride 
comfort criterion without attention to the other criteria and constraints. 

According to eqns. (22)-(24), the system dynamics for xft) i = (2,4) are 
as follows: 


= Si +d jUf +d j+1 u r [i = 2,4) (77) 

Substituting the control inputs (73) and (74) into the open-loop eqn (77) leads to 
the closed-loop dynamics as follows: 

x, + j(d J K i +d j ^k l )x l [l-/i(^.ic, -d j+I A,)] (i-1,2) (78) 

To investigate more specific behaviour of closed-loop dynamics (78), the cheap 
control case i.e. p 7 = p g = 0 is considered. In this case, according to eqns. (75) 
and (76): 


d J K i +d j+l X i =\ 

(79) 

d J K j - d J+1 \ = 0 

(80) 


So, eqn. (78) for the cheap control case takes the following form: 



From the above equation, xft) = 0 for any h > 0. This shows that when there is 
no limitation in using of the applied control forces, the perfect response of zero 
for the sprung mass motion is achieved by ignoring the other criteria 
and constraints. 

4. SIMULATION RESULTS 

In this section, simulation results of the non-linear half-car suspension system 
described in the previous sections are presented to show the effectiveness of the 
proposed control system. Also, the influence of control weighting factors on 
the system responses and control energy are investigated. The physical values of 
vehicle suspension parameters are presented in Table 1. These values are related to 
a typical light commercial car and have been selected near the values of different 
half car models reported in prevalent papers [7,9, 14], 

To simulate the system responses, the road roughness input shown in Figure 2 is 
assumed as several continuous bumps with the following equation: 


x r =0.02(l -sin{5jtt))~ 0.04sin(lOwt), 0.1 s <f <1.3 s (82) 


The road input for the rear suspension is assumed to have r second delay 
according to equation (20). 

To show the adjustability of the controller and the conflictive behaviour of the 
responses, different sets of weighting factors as free parameters are chosen. Firstly, 
the ride comfort criterion, as the primary objective, is investigated by selecting the 
weighting factors as (p 2 = p 4 = 1; p 7 = p g = p x = p 3 = p 5 = p 6 = 0). This case which 
is illustrated in Figure 3, which clearly shows how the accelerations of the front and 
rear suspension stations are perfectly regulated where the sprung mass motion is the 
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Table I. 

Parameters for the case study suspension system 


Sprung mass 

m s 

500 kg 

Front unsprung mass 

m f 

45 kg 

Rear unsprung mass 

m r 

40 kg 

Body moment of inertia 

1 

1200 kgn? 

Distance from the front tire to centroid point 

a 

1.25m 

Distance from the rear tire to centroid point 

b 

1.45m 

Tire stiffness 

k us 

190000 N/m 

Tire damping 

C us 

70 Ns/m 

I st Front spring coefficient 

kfl 

12394 N/m 

2 nd Front spring coefficient 

k fZ 

-73696 N/m 2 

3 rd Front spring coefficient 

kf3 

3170400 N/m 2 

I st Front damper coefficient 

C,| 

1000 Ns/m 

2 nd Front damper coefficient 

C f2 

500 Ns 2 /m 2 

I st rear spring coefficient 

kri 

12150 N/m 

2" d rear spring coefficient 

^2 

-72500 N/m 2 

3 rd rear spring coefficient 

k r3 

3171200 N/m 2 

I st Front damper coefficient 

C rl 

950 Ns/m 

2 nd Front damper coefficient 

C r2 

470 NsV 

Prediction time 

h 

0.1s 



only weighted parameter in the cost function. As a result, the pitch angle magnitude 
is extremely decreased. However, the responses of tire deflections and suspension 
deflections do not converge to steady states in finite time, and they remain severely 
oscillatory. In other words, the road stability of the vehicle is unsafe in this strategy. 
Also, the control forces become large and oscillatory. It is concluded that the 
achievement of best ride comfort is provided at the cost of losing the road holding 
and excessive control forces. Note that, the responses of front and rear suspension 
stations for the uncontrolled system (dotted line shown in Figure 3), are slightly 
different in amplitude and settling time because of the delay and the load transfer 
due to hitting the bumps. 
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Figure 3. Responses of the suspension system when the ride comfort is the only design criterion (p 2 = p 4 = I; 

P 7 = P 8 = °i Pi = Pj = Ps = P 6 = °)- 


As the second case of simulation study, the control inputs are forced to be 
limited in the allowable range by weighting the front and rear control forces in 
the cost function. In this case, a new set of weighting factors (p 2 = p 4 = 1; p 7 = 
p g = I O' 8 ; pj = p 3 = p 5 = p 6 = 0) is considered to limit the control force under the 
maximum value (e.g. 3.5KN). In this case, according to Figure 4, the 
accelerations cannot be regulated exactly with limited control forces and the 
control system cannot achieve the ideal ride comfort. However, the reduction of 
acceleration and also the pitch angle is still acceptable for the ride comfort in 
comparison to the uncontrolled or passive suspension case. The selected 
weighting factors perform the best possible ride comfort in the presence of 
control force constrains. But, the tire deflection responses do not converge to 
steady states in finite time, and remain severely oscillatory as in the previous 
case. In other words, the road stability of the vehicle is still unsafe. This causes 
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Time (sec) Time (sec) 

Figure 4. Responses of the suspension system when the ride comfort with limited control forces is the design criterion 

(p 2 = p 4 ='; p 7 = p 8 = l0_8 ;Pi = p 3 = p 5 = Pt = °)- 


the control inputs to be oscillatory as well. In order to remove the oscillations and 
improve the road holding criterion, the tire deflection and suspension deflection 
of both stations should be weighted in the control law. In this way, a new set of 
weighting factors (p 2 = p 4 = 1; p 7 = p 8 = 5 x 10“ 8 ; p l = p 3 = 5 10“ 7 ; p 5 = p 6 = 

50) was selected. The obtained results in Figure 5 indicate that the suspension 
acceleration and the tire deflection are simultaneously controlled in comparison 
with passive system. Also, the control forces are under the specified maximum 
value (3.5KN). As it is seen, the maximum amplitude of the suspension 
acceleration was reduced to nearly half of the uncontrolled responses in spite of 
applied practical constraints. Therefore, a suitable trade-off between the control 
objectives in terms of ride comfort as primary objective and road holding as 
secondary objective is achieved by considering the restriction on front and rear 
control forces. 
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Figure 5. Responses of the suspension system considering all design criteria (p 2 = p 4 = I; pj = p 8 = 5 X 10 8 ; 
Pi = p 3 = 5xl0- 7 ;p 5 = p 6 = 50). 


To show robust performance of the designed controller, the last set of weighting 
factors has been applied to the controller considering some uncertainties in the 
system parameters. The results of Figure 7 clearly show that the performance of 
the controller did not decrease considerably through uncertainty. In other words, the 
controller has a suitable robustness in the presence of uncertainty. 

The root mean square (RMS) values for all responses from Figures 3 to 6 are also 
presented in Table 2 to offer a quantitative insight into changes in responses. 

Finally, in order to show the results in the frequency domain, the impulse 
response of suspension system is shown in Figure 7. In this case, the road profile is 
modeled as a Dirac delta function and the responses in frequency domain are 
obtained by taking fast Fourier transform (FFT) of each output in the time domain. 
Since the system under investigation is non-linear, applying purely frequency 
domain analysis will not reveal all the system characteristics unlike linear systems. 
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Time (sec) Time (sec) 

Figure 6. Responses of the suspension system in the presence of uncertainties (p 2 = p 4 = I; p ? = p 8 = 5 X 
IO-«; p, = p 3 = 5 x I0- 7 ; p 5 = p 6 = 50). 


However, the results of Figure 7 indicate the improved performance of the 
controlled system in a wide range of frequencies. 

5. CONCLUSION 

In this paper, a predictive control approach has been developed to control a half¬ 
car suspension system model with non-linear elasto-damping forces. Using an 
optimization process, a MIMO non-linear control law has been developed in a 
closed-form which is mathematically tractable and suitable for implementation. 
The results show that a compromise between conflicting objectives in the 
presence of control force limitations can be easily made by regulation of control 
weighting factors. 
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Figure 7. Impulse responses of suspension system in frequency domain (p 2 = p 4 = I; p 7 = p 8 = 5 X 10 8 ; 
p, = p 3 = 5xl0- 7 ;p 5 = p 6 = 50). 


Table 2. 


RMS values of the suspension responses 


(>//4/o X ‘ dt ) 



Front 

suspension 

acceleration 

(*/ =x i} 

Rear 

suspension 

acceleration 

( x ,= x 4 ) 

Front tire 
deflection 
( x i ~*s) 

Pitch 
angle 
( x ,= ) 

Front 

control 

force 

( x , = u) 

Rear 

control 

force 

Uncontrolled 

3.7906 

4.0823 

0.0083 

0.0109 

— 

— 

Controlled in 
Figure 3 

1.0289 

0.9420 

0.0084 

0.0013 

969.9752 

965.9644 

Controlled in 
Figure 4 

2.3090 

2.2036 

0.0076 

0.0041 

473.4736 

503.8095 

Controlled in 
Figure 5 

2.4722 

2.3682 

0.0054 

0.0038 

470.7252 

442.9735 

Controlled in 
Figure 6 with 
uncertainties 

2.4817 

2.3805 

0.0059 

0.0041 

475.2138 

481.3251 
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